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Ballistic persistent currents in disordered metallic rings: Origin of puzzling experimental values 

J. Feilhauer and M. Moskcl 

Institute of Electrical Engineering, Slovak Academy of Sciences, 841 04 Bratislava, Slovakia 

(Dated: December 17, 2012) 

Typical persistent current {hyp) in a normal metal ring with disorder due to random grain boundaries and 
rough edges is calculated microscopically. If disorder is due to the rough edges, a ballistic current Ityp — evp/L 
is found in spite of the diffusive resistance (ex L/l), where vf is the Fermi velocity, I is the mean free path, and 
L ^ I is the ring length. This ballistic current has a simple interpretation: It is due to a single electron that 
moves (almost) in parallel with the rough edges and thus hits them rarely. Our finding agrees with a puzzling 
experimental result Ityp — evp/L, reported by Chandrasekhar et al. [Phys. Rev. Lett. 67, 3578 (1991)] for 
metal rings of length L ~ 1001. If disorder is due to the grain boundaries, our results agree with theoretical result 
Ityp — {evp/ L){1/ L) that holds for the white-noise-like disorder and has been observed in recent experiments. 

PACS numbers: 73.23.-b, 73.23.Ra 



A mesoscopic resistive metal ring pierced by magnetic flux 
($) supports a persistent current IuJ-qII. At zero temperature, 
the ring supports the persistent current / = X]v_e <b h^ 
where /j($) — —dEj{^)/d^ is the single-electron cur- 
rent carried by electron with eigen-energy Ej{^), and Ep is 
the Fermi level fl]. Function /($) is periodic with period 
$0 = h/e, which provides a clear-cut experimental sign of 
the persistent current Uj-Ql]. If the ring is ballistic and pos- 
sesses one conducting channel, the sum ^ Ij changes its sign 
when a new occupied state j is added. Due to the sign cance- 
lation mainly the electron at the Fermi level contributes to the 
sum, and the amplitude of the current is /q — evp/L, where 
vp is the Fermi velocity and L is the ring circumference. If 
the ring is disordered, the size and sign of the cuiTent fluctuate 
from sample to sample due to the disorder fluctuations. It is 
then reasonable to asses a typical cuiTent in a single sample as 
Ityp = (/^)^/^, where (...) means the ensemble average. 

The number of the conducting channels {Nc) in disordered 
metal rings is typically 3> 1 and the rings obey the diffusive 
limit, I <^ L -^ (,, where / is the electron mean free path 
and ^ ~ NJ is the localization length. To estimate Ityp, one 
can assume again that mainly the electron at the Fermi level 
contributes to the sum '^Ij- Since L :^ I, the electron is 
expected to move around the ring by diffusion. Its transit time 
is Tpi = L^/D, where D = vpl/d is the diffusion coefficient 
and d is the dimensionality of the sample. So Ityp ~ e/rp) = 
{1 / d){evp / L){1 / L). A similar result follows from the Green 
functions theory for non-interacting electrons [6], if disorder 
is modeled by a random potential V{r) obeying the white- 
noise condition {V{r)V{r')) oc S{r — r'). One obtains ||6|] 
4^-°^ = 2 X {1.6/d){evp/L){l/L), l^L^ ^. (1) 

Here 2 is the spin factor, d — 1,2, 3, and 1.6 is fromRef. ||7|]. 
The first observation of persistent current in a single metal- 
lic ring was reported f2] for three Au rings of size L ^ 100/. 
The measured currents were ten-to-hundred times larger than 
the result ^; they ranged from ^ O.levp/L to ^ evp/L. 
This huge discrepancy has not been explained yet (see the re- 
views |8, 9]). Other Au rings showed H currents a few times 
larger than result ([!}, and recent measurements of individual 
Au rings ^4^ and Al rings [5] agreed with result dl) very well. 



Why did the similar measurements of diffusive Au lings 



(l2|,|4|] show quite different results, Ityp ~ evp/L and Ityp ~ 
{evp/L){l/L)7 A puzzle |l2[] is why a multichannel disor- 
dered ring of length L ^ I carries the current ^ evp/L, typ- 
ical for a one-channel ballistic ring? This Letter wants to an- 
swer both questions. There is disorder due to polycrystalline 
grains and rough edges H even in pure Au rings. Using a 
scattering-matrix method for non-interacting electrons [Ivl llOll . 
we study typical persistent currents in Au rings with grains 
and rough edges without the white-noise approximation. 

If the disorder is due to the grains, our results agree with 
the white-noise-related formula ([1} and experiments [4, 5]. 
However, if the disorder is due to the rough edges, we find 
the ballistic-like result Ityp ~ evp/L albeit the resistance 
is diffusive (oc L/l) and L ^ I, like in the experiment fj%. 
This ballistic current is due to a single electron that moves 
(almost) in parallel with the rough edges and thus hits them 
rarely. Briefly, result Ityp ~ evp/L in a metal ring of length 
L ^ Hs as normal as the result /typ ~ {evp/L){l/L). Which 
result is observed depends on the nature of disorder 

For simplicity, we study two-dimensional (2D) rings and 
mention the 3D effect at the end. We start with a conductance 
study. Consider a 2D wire [Fig. [T] described by Hamiltonian 

H = -{h'/2m){dl + a^) + C/ [x, y) + V (x, y) , (2) 

where ra is the electron effective mass, U is the grain bound- 
ary potential^and V is the potential due to the edges (see Fig. 
[Hand Refs. ITlllOl]). We connect the wire to two ideal leads - 
clean long wires of width W . In the leads, the wave function 
of the electron with energy E possesses the usual form [12] 

V{x. y) = J::=i [Atix] + A-{x)] sin(Ii^), x < 



^i^,y) = E,"=i mix) + B-{x)] sin(I^), x>L 
where N is the considered number of channels (ideally N 



(3) 



oo), Atix) = a±e±^'="^ S±(x) = &±e 
the wave vector given by equation ^^^ 



±ik„x^ and kn{E) is 
E. The 
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vectors A (0) and B [L) with components A^ 



.N 



(0) and 



B 



n=l....N 



[L), respectively, obey the matrix equation 
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Wire with grain boundaries 



Ax / 



FIG. 1: Model 171. [igi of wires with grain boundaries (a) and rougii 
edges (b). Let d(x) and hi^x) be the j/-coordinates of the edges. 
Then V(x,y) — for d{x) < y < h{x) and V{x,y) — oo 
elsewhere. For smooth edges d{x) — and h{x) — W, other- 
wise d{x) and h{x) fluctuate randomly in intervals (—A, A) and 
{W - A,W + A). The RMS of the fluctuations is 5 = A/y/S. 
The fluctuations occur abruptly with step Aa;. So Ax is the rough- 
ness correlation length. The grain boundaries (a) are modeled as a 
randomly-oriented mutually non-intersecting lines, that is, the angle 
between the given line and the x coordinate is random in a range 
restricted by the presence of the neighboring lines |10]. Each line 
consists of equidistant repulsive dots (plus signs) with potentials 
jS{x ~ Xi)5{y — yi), where (Xi, yi) is the position of the i-th dot. 
Thus [/(a;, 2/) = '^\jilS{x — Xi)5(y — yi). Ifc— > and 7/c is fixed, 
a boundary scatters electrons as a structure-less line-shaped barrier 
independent on the choice of c. If a 2D electron impinges on such 
a barrier perpendicularly with Fermi wave vector kp, it is reflected 
with probability Rq — {'y/c)'^/[kp + (7/c)^], where 7 = rw^jY? . 
The parameters of the grain boundary disorder are Ra (in reality lllll 
Rg ~ 0.1 — 0.8) and the mean inter-boundary distance da- 

where S is the scattering matrix \\T[. Its elements tiE), '>'{E), 
t'{E), and r'{E) are matrices with dimensions N x N. Ma- 
trices t and t' are the transmission amplitudes of the waves 
A"*" and B~ , respectively, and matrices r and r' are the corre- 
sponding reflection amplitudes. In particular, the matrix ele- 
ment tmn{E) is the amplitude of transmission from channel 
n in the left lead into the channel m in the right lead. We 
evaluate S{E) for disorder in figure [T]following Refs. iItIIiOI. 
The wire conductance g (in units 2e^/h) is given by 
the Landauer formula g ^^ ' 



^N. 



Sm^l \tnui{EF)\ 



2krAEF) 



Y.n=iTu, Where T„(Sf) - 
^ ,^ ^ is the transmission probability of 
channel n. We evaluate tmn for many samples and obtain the 
mean transmission (T„) and mean resistance (p) = (l/.g). 

Figure |2] shows our data. The wires with grain boundaries 
exhibit features of the white-noise-like disorder. First, {p) fol- 
lows the usual diffusive law (p) = l/N^ + {2/kFl){L/W), 
shown in a full line. Second, all (T„) are equivalent in the 
sense that (r„) ex 1/L 01311 . Therefore, persistent currents in 
rings made of such wires should agree with equation ([T]l. 

For the wires with rough edges, however, (p) approaches 
the diffusive law (p) = l/N^^f + {2/kFl){L/W) [the fuU 
line in the right panel of Fig. HJa)], where l/N^^^ is the ef- 
fective contact resistance due to the N^^^ open channels. For 
large Nc and small enough Ax we find the universal number 
N^ff ~ 6 < A^c (see [7] and figure SEtb) in [ 14]). These 
open channels dominate also the (r„) data in the right panel 
of figure|2b). Channel n = 1 is almost ballistic ((Ti) ~ 1) 
even for L = 0.2^ ~ 100/, when a few other channels with 
low n show (T„) of the order of 0.1 and all other channels are 
diffusive or localized, with (T„) strongly suppressed OTJ IlSll . 
What are the persistent currents in rings made of such wires? 

Consider the 2D wire in figure [T] but ring-shaped in the 
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FIG. 2: Transport study of the Au wire with grain boundaries [Fig. 
[TJa)] and Au wire with rough edges [Fig. [TJb)]- The parameters of 
Au are m = 9.1 x 10~'^^kg and Ep ~ 5.6eV, other parameters are 
listed. Not to affect the results, TV is usually kept larger than Nc- 
Figure (a) shows the mean resistance (p) versus the wire length L. 
Note that (p) is reduced by the contact resistance 1/Nc and L scaled 
by ^. The localization length ^ is obtained from the numerical data 
for (Ing) (see Refs. (I,[i3]) by means of the fit (Ing) = -L/(_ 
at L 2> C- The full lines show the linear fit of the diffusive regime 
(see the text) from which we obtain the mean free path I. In fact, 
in the right panel one should see four slightly different full lines for 
different Nc ; we show only one of them for simplicity. Figure (b) 
shows (T„) versus L/^ for the parameters indicated by bold arrows. 
For n = 1,2, ... Nc the resulting curves are ordered decreasingly. 

plane of the 2D gas and with the wire ends connected. So we 
have a 2D ring. Since L ^ W, we use the cyclic conditions 






:exp(-z2^$/$o)|f(i,y) 



(5) 



and follow the standard approach H in which the electron 
states of the ring are described by the Hamiltonian of the 
stripe (in our case by Eq. |2]i and by conditions (|5). In other 
words, the ring states are assumed to coincide with the stripe 
states obeying the cyclic conditions (|5]l. Say, in the clean ring 
ip{x, y) ~ sin(^^^)e^*'^''^^, where y is the position along the 
ring radius, cc-axis is bent along the ring circumference, and 
k^{E) = ^(i' + $/$o),wherei/ = 0,±l,±2,.... 

Once the ring states and stripe states coincide (we discuss 
this assumption in detail below), we can estimate the persis- 
tent current in the ring with rough edges intuitively from the 
transmission of the corresponding stripe (the right panel of 
figure |2j)). We do not expect the diffusive result ^ because 
the feature (Ti) ~ 1 at L 3> Hs in contrast with usual diffu- 
sive decay ^ l/L. Assume roughly that (T„) = 1 for n — 1 
and (r„) ^ l/L for all other n. In this model, the chan- 
nel n = 1 contributes by the ballistic current Ityp = evp/L 
while the total contribution from other channels is diffusive, 
Ityp ~ {evF/L){l/L), and negligible for L ^ /. As a result, 
multichannel rings with disorder due to rough edges should 
support at L ^ I typical currents Ityp ~ evp/L. In other 



words, the rough edges scatter all electrons except for a small 
part of those that move (in classical terms) almost in parallel 
with the edges. This small part, composed mainly of elec- 
trons occupying channel n = 1, hits the edges rarely and 
thus moves almost ballistically ll4ll . Eventually mainly the 
electron circulating at the Fermi velocity contributes, as in a 
one-channel ballistic ring. Thus Ityp ~ evp/L albeit L ^ I. 

Now we discuss the assumption that the ring states coin- 
cide for L ^ W with the stripe states obeying the condi- 
tions (|5]l. This standard approach describes the ring by the 
stripe-related Hamiltonian (Eq|2) that ignores the ring cur- 
vature. Consider the clean ring. If the ring curvature is in- 
cluded in the Hamiltonian, it produces the centrifugal force 
which is ex v^ and which pushes the radial wave functions 
towards the outer edge of the ring. Consequently, the radial 
wave functions become localized at the outer edge (especially 
for large ly) and strongly differ from the form sin(^^) even 
for L ^ W. This result is exact in the non-interacting model 
but fails to describe metallic rings, because the localization of 
the radial wave functions at the outer edge produces the in- 
ternal field due to the electron-ion and Hartree-Fock interac- 
tion. This internal field, ignored in the non-interacting model, 
tends to balance the centrifugal force and to delocalize the 
radial waves throughout the ring cross section. Once the bal- 
ance is achieved, the resulting radial wave functions have to 
be close to the stripe-related form sin(^^). Just this is im- 
plicitly assumed in the standard approach that omits from the 
Hamiltonian both the ring curvature and the internal field. In 
reality a small deviation from sin(^^) remains and produces 
the residual internal field balancing the centrifugal force 11611 . 

Consider the standard approach in terms of the semiclas- 
sical paths. Since the x axis is bent along the ring, the x- 
component of any straight-line path in the stripe is bent to fol- 
low the ring curvature; this curvature-mediated orbital effect 
is in fact due to the internal field that balances the centrifugal 
force. The standard approach thus strongly differs from the 
semiclassical-path-based approach [17] that includes the ring 
curvature exactly in the non-interacting model but ignores the 
internal field. As the internal field is ignored, the paths that 
govern the wave functions are the straight lines [17] and the 
radial wave functions are pushed toward the outer edge (this 
is manifested by the straight-line paths that hit solely the outer 
edge [17]). If i ^ W, any straight-line path unavoidably hits 
the ring edges many times [17] , unlike our path that circulates 
almost in parallel with the edges (for further insight see lll4ll ). 

We now verify our estimates of persistent currents by mi- 
croscopic calculations that rely on the standard approach. Us- 
ing Eqs. (|3]l, we write equations (|5]l in the matrix form 

A+(0) 
B-{L) 
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Ring with grain boundaries 



A-(0) 
B+{L) 







Q 







(6) 



where Q is the TV X TV matrix with terms Qq^ = e*^'^*/*0(5Q^. 
Equations (|6]l and Q hold together for discrete energies E — 
Ej{^) which we find for a given ring numerically lllOll . Then 
we find / — -~J2yE <Ef'^^j/'^^ ^^^ ^typ = (P)^^'^, where 
(/^) is averaged over a small energy window at Ep il0lll4ll . 
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FIG. 3: Typical persistent current Ityp in a disordered Au ring versus 
L/l. The ring parameters are shown, <I> — ~-0.25h/e, I has been 
obtained from the wire resistivity (figure |2). The arrows point the 
parameters studied further in Ref. II14II . Symbols are our data, full 

= l.6{evF/L){l/L). 



lines show the formula Ityp"'' 



Figure [3] shows our main results. In the rings with grain 
boundaries, Ityp agrees (at large L) with expected result 
Il^^°^ = l.Q[evF/L){l/L), like in the experiments /3 H/. 
However, in the rings with rough edges, Ityp is systematically 
(not regarding the data fluctuations) close to the ballistic one- 
channel value Iq — evp/L, albeit L ^ I, Nc ^ 1, and 
(p) ex L. All this agrees with the puzzling experiment [A]- 

In the work yO the persistent current ~ /q was observed in 
the Au ring with L ~ lOOZ and W = 90nm. Indeed, the figure 
[3tb) shows Ityp ^ Iq also for L/l ~ 100 and W = 90nm. 
The difference is that in the work H / ~ W^ (^ = 70nm for 
W = 90nm) while our values of I in Fig. [3b) [see also Fig. 
$2a) in [14]] are at least two to three times larger than W; the 
edge roughness alone cannot produce I ~ W. In reality the 
edge roughness coexists with other types of disorder Refer- 
ence |2] did not specify disorder in the measured samples, but 
Webb mentions in Ref. fllSll that the grains in the Au rings of 
work 1 2] were much larger than W (for instance, in Ref. [19|] 
da ~ 8W). The grains with do ^ W are known as bamboo- 
hke grains lH^]. Of course, do :> W and ^ ~ VF |0] means 
I ^ do, which suggests that the grain boundaries were not the 
main source of scattering in the work [^J. If the random grain 
boundaries (or impurities) were the main source of scattering, 
the measured persistent current [2] would be ^ {1/L)Iq rather 
than^ ^0 (c.f. Fig.[3]andRef. [10]). What remains is the edge 
roughness and we have seen that it explains the mysterious co- 
existence of results Ityp ~ /o, L/l ^ 1, and (p) ex L. What 
happens if one adds the bamboo-like grains? 

Since do ^ W, we fit Rq to obtain I ~ W . Figure|4]shows 
such a study for the same W and similar L as in the work |2J]. 
In figure HJa) we see again the diffusive law (p) ex L/l, but 
now I ~ W, like in the work ||2[]. Figure HJb) shows that the 
transmission through channels 1, 2, and a few more is still 
large (between 1 and 0.1), though not as large as in the wire 
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FIG. 4: Transport in Au wires and Au rings with rough edges and 
bamboo-lilce grains. The angle a specifying the orientation of the 
grain boundary is chosen at random from the interval (— ao,ao), 
where qo is the parameter: ao = means the ideal bamboo shape 
with the boundary peipendicular to the wire [20]. The table shows 
all parameters and the resulting I and ^. Figure (a) shows the mean 
resistance (p) as a function of L/l, figure (b) show the transmis- 
sion (Tn) versus L/l for qq — 0. The open symbols in figure (c) 
show the typical current in the ring, Lyp/Io, as a function of L for 
various ao, the full symbols show the corresponding maximum cur- 
rents. Figure (d) shows the Lyp data from figure (c) normalized by 
liyp""^ = 'i--6{evF/L){l/L) and plotted in dependence on L/l. 

with rough edges only [c.f. the right panel of Fig. EJb)]. A 
suppression of the transmission, caused by a combined effect 
of the rough edges and bamboo-like grains, is visible for all 
347 channels. Consequently, I is suppressed as well and we 
have I ~ W. Similarly, the typical currents [Figs, ^c) and 
Ud)] are suppressed compared with the pure edge-roughness 
case [Fig. Ob)], but they still grossly exceed the law I^yp"^ = 
\.&{evF/L){l/L). FigureSfc) shows the maximum currents, 
because the work [2] in fact reported the current amplitudes 
rather than hyp- These amplitudes were between '^ O.l/o and 
~ Iq and roughly the same show our data (full symbols). 

In conclusion, figure |3] naturally explains why the exper- 
iment |12|] shows the result Ityp ~ evp/L and experiments 



confirm the result Ityp ~ {evpJIj){l/L). It suggests 
that disorder in samples of works ||4l |5|] was white-noise- 
like (most likely mainly due to the random grain boundaries), 
while the dominant disorder in [2] was the edge roughness. 

A few remarks at the end. (i) The samples of work [2] were 
3D, but the 3D effects would change our 2D results insignif- 
icantly 11411 . (ii) The step-shaped-roughness model in figure 
[Hb) is universal; our results hold also for any model with a 
smoothly varying roughness 1,1 4.1- (iii) Our results are robust 



against the change of Nc, S, Ax, I, and L for a broad range 
of values. Therefore, the absence of the exact information on 
the nature of disorder in measured samples jM |4l |5| is not cru- 
cial for our conclusions. Anyway, our values of S and Ax are 
realistic (c.f. figures 1 and 2 in Ref. [21]). Experiments that 
would determine Ityp and / in coiTelation with the parameters 
of disorder can be useful, (iv) Note [ 14], that the transmissions 
Tn=i ~ 1 in wires with rough edges have nothing in common 
with the bimodal distribution l/-\/(l — T)T'^, which exists in 
any diffusive conductor [13] and diverges for T = 1. Trans- 
missions T in the bimodal distribution are the eigen-values of 
the t+i matrix, while the meaning of our r„ is different. 
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This supplemental material consists of six sections. In sec- 
tion I we explain in detail why our transport results obtained 
for the roughness model in figure[TJb) hold universally also for 
any other roughness model. In section II, the standard (wave- 
function-based) description of the ring states is compared with 
the semiclassical-paths-based description. In section III we 
provide a further insight into our result Ityp ~ evp/L at 
L/l ::!> 1 by showing supplemental numerical data. In sec- 
tion IV we explain why we do not study the mean current and 
focus solely on the typical current. In section V we explain in 
detail why our 2D study gives the results that hold very well 
also for 3D samples. Finally, in section VI we stress that the 
transmission Tn=i ~ 1 in the wire with rough edges (which is 
responsible for the result Ityp ~ evp/L at L/l ^ 1) should 
not be confused with a well-known general property of any 
diffusive conductor, with the transmission eigen-value T = 1. 

I. On the universality of the roughness model in figurelUb) 

All our transport results for the quasi- ID systems with the 
rough edges, presented in the main text, were obtained for the 
step-shaped-roughness model in figure[TJb). Here we wish to 
point out again that these transport results are universal in the 
sense that they would remain the same also for any model with 
a smoothly varying roughness. 

Evidently, from the technical point of view, the step-shaped 
roughness in figure[TJb) provides a discretization scheme that 
allows to model any smoothly varying roughness by means of 
the very small and very dense steps. Using this approach, our 
calculations from the main text can in principle be repeated 
for any roughness model which is specified by the RMS pa- 
rameter S and roughness correlation length Ax. We expect 
that the obtained transport results will agree with the results 
presented in the main text, if one compares the dependencies 
on the parameter L/^, where ^ is the localization length. This 
expectation is motivated by a few fundamental findings. 

First, the statistical ensemble of the macroscopically- 
identical mesoscopic conductors with impurity disorder is 
known to exhibit the conductance distribution which is essen- 
tially the same (for a given value of parameter i/0 for any 
model of the impurity disorder model; the weaker the disorder 
the better the accord of the conductance distribution of various 
models). Second, it seems that a similar universality (the inde- 
pendence on the specific model of disorder) holds also when 
the impurity disorder is replaced by disorder due to the edge 
roughness. In particular, the conductance calculations in Ref. 
il5l1 . performed for the same step-shaped-roughness model as 
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FIG. Si: (a) The top view on the 2D wire with the rough edges 
generated numerically for two different roughness models. In this 
numerical example the Au wire of width W = 9nm is considered, 
which implies that the number of the conducting channels (Nc) is 34. 
The RMS roughness amplitude S and roughness-correlation length 
Ax are S = 0.87nm and Ax = 0.5nm for the step-shaped rough- 
ness, for the roughness with the Gaussian correlations S — 0.5nm 
and Ax — 1.2nm. Using the approach described in the main text, in 
the former case we obtain the mean free path I — 21nm and local- 
ization length ^ ~ lANcl, and in the latter case we get I = 20.8nm 
and ^ ~ lA9NcL (b) The mean resistance (p) as a function of 
the dimensionless wire length L/^; a comparison for the roughness 
models specified above, (c) The same comparative study as in figure 
(b), but for the mean channel transmissions (Tn); for clarity only the 
data for the first ten conducting channels are presented. 

our model in figure [TJb), gives a quite similar results as the 
conductance calculations in the paper [SMI], performed for 
the smoothly varying roughness with the Gaussian correlation 
function. To demonstrate this universality by means of the 
direct comparison, we have performed the conductance calcu- 
lations for the smoothly-varying roughness with the Gaussian 
correlation function (the model of Ref. [SMI], and we have 
compared them with our results for the step-shaped roughness 
in figurelTJb). 

In figure |2] we show a typical output of such comparative 
study for two Au wires with the same number of the conduct- 
ing channels (Nc = 34), so that one can compare directly the 
individual channel transmission. It can be seen that the in- 



dividual transmissions are in a very good agreement, which 
illustrates the above mentioned universality; note also that the 
individual transmissions for both roughness models coincide 
albeit the values of the parameters S and Ax in the considered 
roughness models are (intentionally) not the same. In addi- 
tion, the universality with respect to the choice of S and Aa; 
within the same roughness model is obvious for all our data in 
the main text and especially from our paper [7]. The last but 
not least, our main result (the ballistic-like persistent current 
Ityp ~ evp/L in figure [3jb)) is universal simply because of 
the absence of the sensitivity to the edge roughness. 

Any user of the scattering-matrix technique can perform 
a similar universality demonstration for other types of the 
smooth roughness. The prize for the use of the smoothly vary- 
ing roughness is a much longer computational time, which is 
crucial especially for the persistent current calculations. Due 
to this reason our study in the main text relies on the step- 
shaped roughness, however, our results are universal as men- 
tioned above. 

II. Comparison of the standard wave-function-based approacli 
with the semiclassical-path-based approach 

The persistent current calculations presented in the main 
text rely on the standard approach. The key assumption of the 
standard approach, justified in the main text, is that the ring 
states coincide for L ^ W with the stripe states obeying the 
conditions (|5]l. Just this assumption allows to describe the ring 
states by the stripe-related Hamiltonian (Eqi2]i that ignores the 
ring curvature. The main text also mentions that the stan- 
dard approach strongly differs from the semiclassical-path- 
based approach |17] that incorporates the ring curvature rig- 
orously (within the non-interacting model) but ignores the in- 
ternal Hartree-Fock-interaction-mediated field balancing the 
centrifugal force. In this approach the semiclassical paths in 
the ring fundamentally differ from those in the stripe, or in 
other words, the ring states and stripe states do not coincide. 
We wish to give a few more comments on the semiclassical- 
path-based approach 1 17]. 

Consider first the clean ring. In spite of the fact that the 
radial wave functions are pushed towards the outer edge of 
the ring by the centrifugal force, the semiclassical-path-based 
approach by Jalabert et al. rtl7ll gives correct result for the 
ballistic persistent current. The problem arises when one con- 
siders the rough edges (see Samohkin's paper lll7ll ). Since 
only the straight-line paths are allowed (neglecting the or- 
bital effect due to the magnetic field), in the annular geometry 
with L ^ W any straight-line path unavoidably hits the ring 
edges many times when it makes one trip around the ring | ITJ. 
On the contrary, the stripe geometry allows also very long 
straight-line paths that are (almost) parallel with the stripe 
edges and therefore do not feel the edge roughness. Within 
the approach of Refs. [17], these very long ballistic paths are 
changed in the annular geometry on the paths that hit exclu- 
sively the outer edge of the ring and become scattered by the 
edge roughness. On the other hand, in our standard approach 



these paths remain parallel with the ring edges and carry the 
ballistic-like current. This fundamental difference is due to 
the fact that the standard approach implicitly incorporates the 
balance between the centrifugal force and internal Hartree- 
Fock field (see the main text), while in the semiclassical-path- 
based approach 11711 the centrifugal force remains unbalanced. 

We also wish to address another fundamental difference. 
Samokhin fl7] assumed that the straight-line path that hits 
the ring edge is reflected diffusively no matter what is the 
incidence angle (the angle between the path and the edge). 
Specifically, the probability of the diffusive reflection from 
the edge is described by the Fuchs coefficient which is equal 
to unity for all incidence angles. It should be mentioned that 
a realistic probability of the diffusive reflection, derived by 
Soffer and Ziman [SM2] for a free wave impinging the sur- 
face with uncorrected roughness, strongly depends on the in- 
cidence angle. In particular, it is equal to unity for perpendic- 
ular incidence but approaches zero for small incidence angles. 
Samokhin lll7ll found in the ring with rough edges the diffu- 
sive persistent current. However, he would certainly not find 
such diffusive current, if the Fuchs coefficient is replaced by 
the angle-dependent probability of the diffusive reflection due 
to Ziman and Soffer [SM2]: Owing to the (almost) specular 
reflections at small angles, his result would most likely be- 
come more similar to our ballistic prediction. In our paper the 
correct angle dependence of the edge roughness scattering is 
included microscopically in the scattering matrix method. 

Indeed, the tendency to a specular reflection at small angles 
is manifested by the channel transmission r„. Let us look 
at the right panel of figure 2b in detail. Semiclassically, the 
channel number n corresponds to the angle between the semi- 
classical trajectory and edge, and n = I corresponds to the 
smallest nonzero semiclassical angle allowed by the quantum 
confinement. Consider, say L ~ 0.25^ ~ 1201. In case of 
the diffusive reflection, for L/l = 120 one should observe 
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FIG. S 2: (a) The mean free path / and (b) effective number of the 
open channels N^ in the wire with rough edges, both plotted in 
dependence on the roughness correlation length Aa; for the param- 
eters as indicated. These data were extracted from the data for (p) 
versus L by means of the fit (p) = l/N^^^ + {2/kFl){L/W) as is 
explained in Fig. |2ja) and in the main text. For simplicity, 5/W is 
kept nearly the same (~ 1/10) for each set of 5 and W . 



T„ ~' l/L ^ 1/120 already for n — \. However, this is 
not the case; the right panel of figure 2b shows that r„ is be- 
tween 1 and 0.1 for n = 1, 2, . . . , 6. In semiclassical termi- 
nology, transmissions for "angles" n — 1, 2, . . . , 6 exceed the 
diffusive value ~ l/L ~ 1/120 by one to two orders of mag- 
nitude. This means that the motion within these channels is 
much more ballistic than diffusive, although (semiclassically 
speaking) the electrons in these channels experience multi- 
ple collisions with the edges. The diffusive values of Tn are 
reached only if n is large enough, that means for the large 
enough semiclassical angles. 



III. Supplemental numerical data 

The figure S|2] shows the electron mean free path / and the 
effective number of the open channels N^f^ in the wires with 
rough edges, extracted from the linear fit of the numerical data 
for (p) versus L (see the discussion of figure |2] in the main 
text). Two features in the figure are worth noticing. 

First, the I versus Aa; dependence in figure S|2a) demon- 
strates clearly that the minimum mean free path due to the 
edge roughness scattering is always a few times larger than the 
wire width W . In other words, the edge roughness alone can- 
not explain the observation I ~ W, reported in the experiment 
yO- To explain the observed persistent currents / ~ evp/L 
jointly with observation / ~ w, the edge roughness has to be 
combined with the bamboo-like grains, as is demonstrated in 
figure|4]in the main text. 

Second, the N^^^ versus Aa; dependence in figure SHIb) 
demonstrates that N^^^ is a universal (A^c-independent) num- 
ber of the order of 10 if Aa; is small enough and Nc large. The 
universal N^^^ has been discovered in Ref. |7|], here the uni- 
versality is demonstrated for Nc as large as 347. Of course, 
the rings made of such wires have to possess the same N^^^ ■ 
Now we demonstrate that this is indeed the case. 

In figure S[3]we show how the typical current in the ring 
with rough edges depends on the number of channels (N) 
considered in the simulation. It is (roughly) A^-independent 
for N > 10, no matter how large Nc is. In other words, 
the currents ^ Iq in rings with rough edges are almost ex- 
clusively carried by the open channels n — 1,2,..., Nc^^, 
where Nc^^ ^ 10 for any value of Nc- 

Moreover, as the transmission of these open channels is 
large for L/l 3> 1 [in particular, (Ti) ~ 1, as is shown 
in the right panel of figure |2b)], one also understands why 
the ring with rough edges supports the ballistic-like current 
Ityp ^ Iq rather than the diffusive current Ityp ^ {1/L)Iq. 
The current Ityp ^ /q exists because the rough edges scat- 
ter all electrons except for a small part of those that move 
(classically speaking) almost in parallel with the edges. This 
small part, composed mainly of electrons occupying channel 
n = 1, hits the rough edges rarely and therefore moves al- 
most ballistically. The resulting current, Ityp ^ evp/L, is 
due the electron that circulates with Fermi velocity, because 
contributions from other electrons tend to cancel like in a true 
single-channel ballistic ring. 
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FIG. S 3: Typical persistent current Ityp in the ring with rough edges 
as a function of the total number of channels (TV) considered in the 
simulation. The same parameters and symbols are used as in figure 
|3jb), the considered ring lengths are shown as L/l. 



Since (Ti) ^ 1, one could naively think that the value 
Ityp ^ Iq will survive also if one chooses N as small as 
N — 1. Figure $3] shows that this is not the case. For in- 
stance, in the ring with Nc — 347 and L/l = 120 the current 
approaches zero just for N —>■ 1. This is easy to understand: 
Once the channel n = 1 cannot communicate with other chan- 
nels, the transmission (Ti ) ^ 1 tends to be suppressed to zero 
by Anderson localization, present in any sufficiently long ID 
disordered system. Communication with a few other channels 
is needed to restore (Ti ) ^ 1 and to obtain Ityp ^ Iq . 

To provide further insight, figure S2] shows the sample- 
specific currents in two selected rings from figure |3] (bold 
arrows) and in a clean ring. Figure SSJa) shows the de- 
pendence Ij versus Ej, figure ^Ifb) shows the total current 



E 



/, versus Ep. Evidently, the ring with rough 



edges exhibits remarkably larger cuiTents than the ling with 
grain boundaries, albeit both lings are of the same size and 
posses the same value of I. 

Figures SUc) and $5d) focus on a small energy window 
below the Au Fermi level. One can see that Ij in the ring with 
rough edges exhibits sharp peaks with the sign alternating and 
oscillating with period AE = 2TrhvF/L. This period is twice 
the inter-level distance in the ballistic single-channel ring, 
which suggests that the peaks are due to the quasi-ballistic 
channel n = 1. [We recall that (Ti) -- 1 also for L/l > 1, 
as is shown in the light panel of figure |2lb).] However, the 
height of the peaks is affected also by other channels, because, 
as discussed above, channel 1 cannot keep (Ti) ^ 1 without 
communicating with a few other channels. 

In figure S|4jd) one can see that in the ring with rough edges 
also the total current I{Ep) oscillates with period AE. The 
amplitudes of the total current are close to Iq, and therefore 
the typical currents of size ^ Iq appear in figure|3jb). 

In fact, aheady the data for the clean ring show I{Ep) os- 
cillating with period AE. However, the amplitude of / is 
^ a//Vc2/o and the amplitude of /„ is 2/o, where the factor of 
2 is due to the spin. Evidently, the rough edges reduce / from 
^ -\//Vc2/o to '^ Jq, but they do not change the oscillation 
period set by the clean ring. Note that also the ring with grain 
boundaries exhibits the oscillating persistent current. These 
oscillations are chaotic and correlated with correlation length 
~ {l/L)AE, predicted |6] for the white-noise-like disorder. 
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FIG. S 4: Persistent currents in a ring with grain boundaries, a ring 
with rough edges, and a clean ring for the parameters marked by the 
arrows in Fig.|3] for L — 375 nm, and for $ = — 0.25/i/e. For both 
disordered rings, the considered parameters ensure 1{Ef) ~ 21 nm 
at the Au Fermi level (Ep ~ 5.6 eV). Figure (a) shows the single- 
electron current Ij versus the eigen-energy Ej . Figure (b) shows the 
total current I = X^y^ <g Ij obtained by summing the currents in 
the figure (a) for Ep varied from to 5.6 eV. Figures (c) and (d) show 
the same data as the figures (a) and (b), but for a small energy window 
below the Au Fermi level. The data are scaled by /o — evp/L, the 
data points are connected by full lines which serve as a guide for the 
eye, the bars depict the energy increment A_E = 2TThvF /L. Figure 
(e) shows the typical current Ityp = (1^)^". Averaging over the 
energy window in figure (d) gives the values shown by dashed lines: 
hyp/Io — 1.6(//L) for the ring with grain boundaries, hyp/ lo — 
0.5 for the ring with rough edges, and Ityp/Io — \fNc for the clean 
ring. The circles show the data obtained by varying the number of 
channels, A^, from A^ = 1 to A^ > A^c (here Nc = 34). 

Figure S|4je) shows the typical current. The dashed lines 
show the values of Ityp obtained from the data in figure SlUd), 
the circles show Ityp in dependence on N . For all three 
rings one sees, that the circles approach with raising N the 
iV-independent value (the large N limit) represented by the 



dashed line. It can be seen that a reliable estimate of Ityp in 
the ring with grain boundaries requires N > Nc, while for the 
ring with rough edges one only needs N ^ 10 no matter how 
large N^ is. This is due to the effective number N^^^ ~ 10, 
as has akeady been explained in the beginning of this section. 



IV. The problem of the mean persistent current 

In our present work we have focused on the typical current 
and we have not discussed the mean current. The sign and 
amplitude of the mean c_urrent, measured in the pioneer ex- 
periment by Levy et al. JjJ] is another puzzling problem in the 
field. As can be seen from our description of the scattering 
matrix method, we can in principle provide also numerical 
data for the average persistent current. There are however a 
few serious reasons why our present manuscript is not focused 
on the average current. 

First, the problem of the average current has been addressed 
by Bary-Soroker et al. [SM3] who attempted to explain it 
within the interacting electron model. However, these authors 
did not address the problem of the typical current, and the 
experiments 141 15|] showed, that the typical current is most 
likely not affected by electron-electron interaction and should 
be tractable within the non-interacting model. Just these rea- 
sons lead us to focus on the problem of the giant typical cur- 
rent and to solve it within the single-electron model. 

Second, a complete scattering-matrix study of the ampli- 
tude and sign of the average current would require (perhaps) 
ten to hundred times more computational time than the study 
of the typical current, presented in our present manuscript. We 
have therefore decided to focus on the problem of the typical 
current, and akeady this problem was computationally cost. 



V. On the robustness of our 2D results against the 3D effects 

All our transport data in the main text were obtained within 
the 2D model depicted in figure [T] while the experimental 
samples of reference [2] were three-dimensional. Here we 
want to point out in detail that the extension of our 2D study 
to 3D (replacement of the rough edges by rough side-walls) 
would not change our major results remarkably. The effect 
of 3D on our 2D results can be estimated easy without any 
explicit calculation as follows. 

In our 2D wire ([r|(b)) only the edge roughness scattering is 
considered, while in the 3D wire of reference [2J the rough- 
ness scattering is in general due to the wire edges (side walls) 
and due to the top and bottom surfaces in addition. In spite of 
this difference it is evident that the 3D sample still preserves 
the key feature of our 2D model: The electrons which occupy 
the ground ID channel (now the channel with quantum num- 
bers Uy = 1 and Tiz — 1, where z is the vertical direction) still 
move almost in parallel with the sample edges and sample 
surfaces, and therefore avoid the roughness scattering quite 
similarly as in the 2D case. Due to this key feature, the almost 



perfect transmission of the Fermi electron in the ground ID 
channel, found for the 2D wire (Fig. 2), has to persist also in 
the 3D wire. Owing to this ballistic Fermi electron, the 3D 
ring has to carry the persistent current Ityp ~ evp/L also for 
L/l >> 1, just like the 2D rings in figure|3lb). 



Furthermore, one can easy see that the roughness-mediated 
scattering in 3D does not modifies the mean free path I in com- 
parison with 2D remarkably. As akeady mentioned, in the 3D 
wires one should consider also the roughness-mediated scat- 
tering from the top surface and bottom surface. However, in 
real normal-metal 3D wires the roughness amplitude (RMS) 
of the top and bottom surfaces is usually of the order of one 
lattice constant (^ 0.5nm; see e.g. the paper [SM4]), which is 
one order of magnitude less than the roughness amplitude at 
the edges (RMS ^ 5nm - lOnm; see the experiment of 1,21,1 , 
and our manuscript). Since the roughness-limited mean free 
path is proportional to the square of the RMS (see e.g. our 
paper |7] or references cited therein), the effect of the top and 
bottom surfaces on the mean free path is two orders of magni- 
tude weaker than the effect of the edges. In other words, it is 
very likely, that in the 3D wires of reference [2J the roughness- 
mediated scattering is mainly due to the wire edges, while the 
reflections from the almost smooth top/bottom surface are (al- 
most) specular and thus do not affect the electron transport. 
The case when the roughness of the top and bottom surfaces 
is comparable with the edge roughness is experimentally un- 
like: The roughness of the edges is due to the limitations of the 
electron lithography and lift-off, due to the roughness of the 
resist side walls, etc. (see [21]), while the origin of the surface 
roughness is quite different. However, even in this unlike case 
the rough surfaces would marginally affect the quasi-ballistic 
channel Uy = I and n^ = 1 and they would only reduce the 
edge-roughness-limited mean free path by a factor of ^ 2, as 
can be estimated from the Mathiessen rule. 



Finally, unlike our 2D wire in figure [T] (b), the edges of 
the 3D wire are in fact the side walls and the edge rough- 
ness at such side walls in general scatters the electrons also 
into the vertical (z) direction, in addition to the in-plane scat- 
tering considered in our 2D model. This may decrease the 
roughness-limited mean free path say by a few tens of per- 
cent, but cannot affect the fundamental feature (the ballistic- 
like motion in the ground ID channel) responsible for the per- 
sistent current Ityp ~ evp/L also for L/l >> 1. 



In principle, our scattering-matrix approach allows to pro- 
duce also the numerical data for the 3D samples (we already 
did so for the wires with the grain boundaries in reference 
OlOD), albeit not for the samples as large as those used in the 
experiment ^. However, the computational time would be 
quite huge and we do not expect any new features in compar- 
ison with our 2D results. 



VI. Feature T„=i ~ 1 in the wire witli rough edges and T = 1 as 
a general feature of any diffusive wire: Two different things 

We recall that the ring with the rough edges exhibits the 
persistent current Ityp ~ evp/L at L ^ I (figure |3lb)), be- 
cause the constituting rough wire exhibits the transmission 
T„=i ~ 1 (right panel of figure |2b)). We want to stress that 
the transmission Tn=i ~ 1 in the wire with rough edges has 
nothing in common with the well-known bimodal distribution 



l/-\/(l — T)r2, which exist in any diffusive conductor lll3ll 
and diverges for T = 1. Note that the transmissions T in 
the bimodal distribution are the eigen-values of the t^t ma- 
trix ll3ll . while we speak about T„ = X]m=i l^n.mP' which 
are the diagonal elements of the i+t matrix. In other words, 
the channels corresponding to the eigen-values T in the distri- 
bution l/-y/(l — r)T2 are the eigen-states of the f+t matrix, 
while the channels corresponding to our diagonal elements T„ 
are the usual plane-wave states. To avoid misunderstanding by 
readers who are accustomed to the latter definition, we add a 
few more remarks. 

The bimodal distribution l/-\/(l — T)T'^ as a general prop- 
erty of any diffusive conductor with white-noise-like disor- 
der 1I13I1 coexists with the diffusive persistent current Ityp ~ 
{evF/L)(l/L) in the corresponding disordered ring [J61. In 
other words, the eigenvalues T = 1 in the bimodal distribution 
do not cause any ballistic persistent current. The reason why 
the current is diffusive in spite of T = 1, is most likely that 
the eigenvalue T = 1 does not necessarily mean the ballis- 
tic motion (a well known example is the perfect transmission 
in case of resonant tunneling). The situation becomes funda- 
mentally different for disorder due to the rough edges. In this 
case the eigen-values Tn still follow the bimodal distribution 
l/-\/(l — T)T'^, however, this has nothing in common with 
the ballistic-like persistent current found by us. The ballistic- 
like current is due to the appearance of the diagonal element 
Tn=i — 1: Namely, any wire from the statistical ensem- 
ble of wires with rough edges exhibits the diagonal element 
Tn=i — 1 independently on the choice of the Fermi energy 
and wire length. These features are the attributes of ballistic 
electron motion within channel n = 1. Indeed, it is easy to 
check in our simulation, that the electron plane wave entering 
the wire in the channel n = I remains (almost) unscattered 
between any two successive scatterers inside the disordered 
region. As a result, the ring made of such wire supports the 
persistent current with typical size dominated by the ballistic 
channel n = 1, that is, Ityp ~ evp/L. In summary, the reason 
for appearance of /typ ~ eui?/Lis the ballistic behavior of the 
diagonal element T„=i; the fact that the bimodal distribution 
gives eigenvalues T = 1 is irrelevant. 
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